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ABSTRACT 


In this study we describe an innovative method to determine potential sites 
for optical and infrared astronomical observations in the Andes region of northern 
South America. The method computes the Clear sky fraction (CSF) from Geosta¬ 
tionary Observational Environmental Satellite (GOES) data for the years 2008- 
12 through a comparison with temperatures obtained from long-term records of 
weather stations and atmospheric temperature prohles from radiosonde. Griteria 
for sky clearance were established for two infrared GOES channels in order to 
determine potential sites in the Andes region of northern South-America. The 
method was validated using the reported observed hours at Observatorio Nacional 
de Llano del Hato in Venezuela. Separate GSF percentages were computed for 
dry and rainy seasons for both, photometric and spectroscopic night qualities. 
Twelve sites with Eve year averages of GSF for spectroscopic nights larger than 
30% during the dry seasons were found to be suitable for astronomical observa¬ 
tions. The best site with (220±42) spectroscopic clear nights per year is located 
in the Andes of Venezuela (70°28’48”W, 9°5’60”N) at an altitude of 3480 meters. 
Lower quality regions were found in Sierra Nevada de Santamarta and Serrama 
del Perijd with (126±34) and (111T27) clear nights per year, respectively. Sites 
over the Andes are identified in Norte de Santander with (107±23) and in the 
north-east part of Boyacd with a mean of (94±13) clear nights per year. Two 
sites at low latitude located in Ecuador with more than 100 clear nights per 
year and with similar seasonal GSF percentages were also identihed. Five year 
evolution suggest a possible correlation between the lowest percentages observed 
during the rainy seasons of 2010 and 2011 with positive values of the Southern 
Oscillation Index. 
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Subject headings: Astronomical Instrumentation. 



1. Introduction 


The use of images of the earth in different wavelengths from infrared cameras on board 
of geostationary satellites is a powerful tool for the identihcation of new potential sites 
for astronomical observations in optical, infrared and submillimeter wavelengths (Erasmus 
& Sarazin 2002; Cavazzani et al. 2012; Cavazzani & Zitelli 2013). Sites in remote areas 
with complex relief and scarse or null meteorological records can be identihed by previous 
computation of the Clear Sky Fraction (CSF). As the cloud covering is a function of the 
height of the atmosphere, images from satellite are useful due to the detectors on board 
measure the emissivity of the atmosphere at different layers. 

Several studies for the determination of the CSF are found in the literature. Soden 
& Bretherton (1993) conducted an analysis of upper tropospheric humidity and clouds 
using GOES observations. Ending a correlation between humidity and cloud cover at high 
altitudes. Erasmus & Sarazin (2002) were among the first to sucesfully use long-term 
satellite data for CSF computations for a wide range of altitudes. Erasmus & van Rooyen 
(2006) and Della Valle et al. (2010) reproduced with high precision the CSF records of 
ground data at la Palma with a model based on satellite data. Cavazzani et al. (2011,2012) 
found percentages of clear nights at five astronomical observatories around the world using 
long-term GOES data. Hidayat et al. (2012) used 15 years of GOES infrared images in 6.7 
and 10.7/im bands to compute the percentage of clear sky over a wide region in Indonesia. 
Recently, Marin et al. (2015) computed the CSF at high altitude sites by using a new 
clearance algorithm based on global meteorological models. CSF values were then used for 
precipitable water vapour estimations at two astronomical sites in the north of Chile. 

Even though most of the methods are based on the comparison of radiances, fluxes or 
temperatures at each pixel of the satellite image with threshold values previously defined 
there is no standard methodology for CSF computations. Cavazzani & Zitelli (2013) used 
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monthly thresholds dehned as the maximum radiance at clear ni^ts along the month. 
Thresholds may also be dehned by using an atmospheric model (j9|). In this study, we 
introduce a novel methodology involving the use of monthly thresholds of temperature at 
different altitudes from the surface up to the upper troposphere. Our thresholds dehnitions 
were based on the altitude of the terrain, long-term records of weather stations and 
radiosonde data archives. 


The thresholds dehnitions are based on monthly mean records of both, air temperature 
at the site and vertical temperature prohle between 8 and 10 km of altitude which make 
more precise the predictions for clear covering than those involving atmospheric models 
for thresholds computations. Our methodology was applied in northern South America 
over the tropical region covering Colombia, western Venezuela and northern Ecuador. Even 
though the zone is far from being interesting for infrared and submilimeter astronomy due 
to the requirement of excellent weather conditions, many universities and institutes from 
those countries are planning to develope optical observatories. The countries are located 
within the intertropical convergence zone (ITCZ) showing complex topography and with a 
climate governed by a bimodal annual cycle which results of the double passage of the ITCZ 
for the region (Poveda et ah 2011). The Andes affect the climate by blocking zonalflows, 
influencing regional wind patterns and precipitation and hosting diverse micro-climates. 
The main goal of this study is to seek prospective locations with middle quality CSF values 
over the Andes of those countries. 

The paper is organized as follows. In Section 2 we describe the satellite data used in 
this study. Threshold temperatures calculations are presented in Sections 3 and 4. The 
methodology and validation are discussed in Sections 5 and 6 while in Sections 7 and 8 we 
present the results and the selected sites. 
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2. Satellite based data 


The Geostationary Operational Enviromental Satellites GOES 12 and 13 are eqniped 
with instrnmentation capable of captnre images in several infrared bands. For onr pnrposes 
we concentrate onr attention on B3 and B4 bands centered at 6.7 and 10.7/rm, respectively. 
Images in B4 are strongly affected by emission and absorption of water vaponr fonnd at 
middle altitndes (~4 km) whereas layers of the npper trophosphere between 8 and 10 km 
are sensed by images in the B3 band. Other bands as those centered in 13.3/im provide 
information abont fog and very low atmosphere characteristics. 


Since the data extends for decades these satellites are commonly nsed for long-term 
stndies. Temporal domaine for this stndy begins on the 1th of Jannary 2008 and ends on 
the 31st of December 2012 with a few missing data lower than 7% in the worst cases. GOES 
missing data has different canses: technical tronbles in GOES-12 (from December 15th 2008 
to Jannary 6th 2009); when GOES-13 replaced GOES-12 (from April 14th to 28th 2010); 
and due to anomalies in the operational system (from September 23th to October 18th 
2012). We assumed the clear sky fraction as null for the days with GOES missing data. 


GOES 12-13 images are freely available through an electronic library of NOAA^J in 
format GVAR (GOES Variable Format). We selected a region of the image to download 
covering a big part of Golombia, the western part of Venezuela and the north of Ecuador. 
For each night we used available observed time at the three different hours 00:15:15, 
06:45:13 and 09:45:14 UTG in order to cover the entire night. GVAR is an instrumental 
radiance given in counts. We used the library of idl routines provided by the National 
Glimatic Data Genter (NGDG) for converting the counts to radiance at each pixel. Usual 
calibration of the images was done using the scaling bias and gain values available from 


^National Oceanic and Amospheric Administration 
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NCDC websitqj in order to get the radiance in mW/[m‘^srcm 


-n 


Brightness temperatures at the instants i = 1,2,3 of the night are computed 
from calibrated radiance values by assuming black body emission. Conversion to actual 
temperatures is done using the relation Tj = a + where a and (3 are conversion 
coefficients dependent on each band (Weinreb et ah 1997, Rossow et ah 1985). Reported 
values available from NCDC website, for GOES 12 satellite are [3 =1.01 and a = —4.76K 
for the B3 band and (3 =1.00 and a = —0.36K for the B4 band. The differences between 
the actual and brightness values increase with decreasing temperature but usually are lower 
than 4%. For instance, a brightness temperature of lOOOOK corresponds to temperatures 
of 1095K and 9999.6K for the B3 and B4 bands, respectively. For most of cases those 
differences are negligible for the CSF computations. We indicate actual temperatures in B3 
and B4 bands as and TB 4 ^j, respectively. 


The spatial resolution of GOES images is 4x4 km in both bands and thus a single pixel 
subtends only a small held of view from the ground compared with the total held of view 
available for an observer at ground. As the validation of our method involves comparison 
with log-books of the Observatorio Nacional de Llano del Hato in Venezuela (hereafter, 
CID) we degraded the original resolution in order to make more real comparisons. Other 
authors have used similar or even lower resolutions. Cavazzani and Zitelli (2012) computed 
GOES temperatures for each site over a 1° x 1° sub-matrix which corresponds to a linear 
projection of about 100x100 km centered on the coordinates observatories at ground. We 
used a Delanuay interpolation of actual temperatures Tes* and TBA,i on a grid of 51x51 
pixels obtaining a resolution of 28x18 km which we consider wide enough for validation 
with ground data. 


^http://www.ncdc.noaa.gov/ 





3. Threshold temperature for middle-high clouds detection 


Detection of mid-altitude clouds at a height of 4000 meters above a pixel location of the 
image in B4 band and at a hxed instant during the night requires the previous evaluation 
of a threshold or reference temperature. In this study we compute monthly thresholds 
of temperature based on long-term records of weather stations located inside the region 
covered by the GOES B4 and B3 images. 


We used long-term records of meteorological data obtained from the Instituto de 
Hidrologia, Meteorologia y Estudios Ambientales, hereafter IDEAM. The IDEAM data 
archive comprise periods from 7 to 40 years of monthly mean values of temperature and 
precipitation for a set of weather stations in Colombia. The complete archive contains 
400 records at altitudes between 400 and 4700 meters above sea level covering 25% of the 
country. The archive presents a high number of gaps, especially those of temperature. We 
considered only long-term records of temperature with a maximum of one monthly absence 
obtaining the reduced number of 63 weather stations irregularly distributed in a region 
comprising latitudes between 3° and 7° North, and longitudes from -76° up to -72° West. 
We will refer to this region as the control region defined by the spatial distribution of the 
IDEAM long-term historical records. The longest records of these selected stations have 
more than 40 yr as is shown in the right panel of Figured] Locations are widely distributed 
in different altitudes with a high number of them close around Bogota (BOG) (h=2600 
meters) as shown in the left panel of the same figure. 


The average of monthly mean records at each station for periods spanning at least 
7 years is expected to be normally distributed around a central value < Tg > being very 
likely that the monthly mean values are within Sag of the distribution, where ag is the 
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standard deviation. The coldest temperatnre at each station may be then dehned as 
Tmin =< Tg > —3(Js. However, an additional correction term mnst be applied since the 
measnrements do not represent valnes during the night. We recall that the IDEAM archive 
consists only of mean temperatures thus we must correct them for nightime observations. 
The monthly threshold temperature at each station may therefore computed as : 


Tq — 2 [^min 


( 1 ) 


where Tmin,tr is the trend of the minima records of temperature given by short-term records 
of weather stations indicated with a solid line in Figure |2] and represents the correction 
for night-time observations. We note that these thresholds correlate very well with the 
altitude of the station as illustrated in the same figure suggesting that monthly threshold 
temperatures may be completely determined from the altitude of the stations. From linear 
hts for each month we established monthly relationships between threshold and altitude 
of the form: To(h) = A x h + B where A and B are constants for each month with values 
indicated in Table [H Slight differences along the year are observed for these two parameters 
with slope or lapse rates variations of the order of 10“^ °C and negligible variations in the 
intercepts. 


3.1. Kriging interpolation of altitudes 

Thresholds in the control region at locations where no weather stations are available 
require interpolation of the altitudes. We used the geostatistical method called Kriging, 
which is a spatial interpolation method that in its simplest form allows to compute altitudes 
at locations with no data, based on the spatial autocorrelation. This means that two close 
stations have similar values of altitude whilst two stations farther apart do not. Briefly, the 
height at the location sq refered as h{so) is given by a weighted combination of the IDEAM 


all records of altitude : 


400 

h{so) = '^Wi{so) ■ h{si) = Xq ■ h (2) 

i=l 

where Aq is the vector of Kriging weights Wi and h is the vector of available samples 
of altitudes. Our sample is dehned by the altitudes of all available records of weather 
stations. The system of equations ([2]) is equivalent to the matrix form Aq = C~^ ■ cq where 
C represents the matrix of covariances between all sampled values and Cq the vector of 
covariances between sampled points and the point in which the value is assessed. The 
vector of covariance is obtained from the semivariogram 'y{d) dehned as: 

7(d) = T|(;4si) - h{s, + d))^] (3) 

where d represents the unit distance or lag from which semivariance is determined. 
The two semivariograms indicated in Figure [3] were computed for lags up to 200 km along 
two mutual perpendicular directions. Both semivariograms present similar spatial variances 
and are well described by a spherical model (dashed line) with a nugget of 0.048, sill of 
1.01 and range of 120. The spherical model was found as the most suitable for modeling 
the covariances ([3]) and for solving the weights in the system of equations (|2]) by using 
GSL-GNU numerical libraries for matrix inversion. A comparison of the interpolated 
heights with real altitudes from all-term IDEAM weather stations is shown in the Figure 
m Gross validation of the Kriging interpolator was done by hiding each IDEAM data point 
during the computation of the weights and then comparing the prediction with the real 
value. Differences with an ideal Kriging interpolation (solid line) are within 3(j which we 
consider enough good for our purposes. 


Based on Kriging results for altitude interpolation in the control region, we computed 
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the monthly thresholds using the values appearing in the Table [H Results are very similar 
for all the months with small differences due to the bimodal annual cycle. Threshold 
temperatures for January with resolution of 28x18 km are shown in the panel b) of Figure 
m The coldest region (0-2)°C corresponds to the Nevada del Huila the highest snow peak 
in the central branch of the Andes. Thresholds between (2-7)°C are observed in the eastern 
Andes comprising Villa de Leyva (VDL), Pisba (PIS) and Macaravita (MAC). Warm 
regions with thresholds above 20°C characterize the valley of Magdalena river between the 
central and eastern branches of the Andes and the southestern regions on zero altitude 
zones. 

For comparison purposes, we show in panel a) the thresholds obtained using Delaunay 
triangulation of the 63 IDEAM long-term monthly records. We note that the spatial 
domaine is constrained in this case and defined by the triangulation of the 63 irregularly 
distributed long-term IDEAM weather stations. Differences between both interpolations 
are signihcant mostly in regions surrounded by strong topographic gradients in which 
a better estimation of the threshold temperatures using Kriging is expected. This is 
particularly noticed in the poor IDEAM sampled region centered on PIS where a maximum 
difference of ±8°C is observed. Other regions with low number of weather stations are the 
south of Bucaramanga (BUG) and the south-east of Medellin (MED). On these regions 
the combination of the lack of weather stations and strong gradients in altitude make 
the Delanuay interpolation unrealistic. Mean differences of ±3°C over the eastern Andes 
are atribuited mainly to strong gradients as the region is well sampled by the IDEAM 
data archive. Thresholds obtained from altitude are more precise than those obtained by 
interpolation of weather station records and are used as a reference for cloud detection in 
the B4 band. 


Uncertainties in the threshold temperature determination are computed using the 
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sigma error of the Kriging interpolation. The largest errors up to 6°C are observed over the 
sea level locations and into the valleys while errors over the Andes are always below 1.5°C 
as shown in panel c) of Figure 01 For the control region the thresholds are computed with 
an uncertainty of ±2°C. 

With the aim of computing monthly threshold temperatures over a wider region 
including the Observatorio Nacional de Llano del Hato where log-books of cloud cover 
fraction are available, we apply monthly threshold-altitude relations over what we denoted 
as extended region covering latitudes between 0 and 13N and longitudes between 70 
and 78W. A set of 3000 altitudes over this region were obtained from Google-Maps and 
interpolated using the methodology discussed above. Monthly threshold temperatures 
between 21 and 24°C are observed across Caribbean coastal regions and in the valleys of 
the rivers Magdalena and Cauca in panels of Figure [5l In the proximity to the Andes, 
thresholds get colder specially in the highest mountains indicated with purple in the hgure 
such as the four highests locations over the Colombian and Venezuelan Andes highlands 
with threshold temperatures below 0°C. Uncertainties in threshold temperatures over the 
extended region are assumed to be equal to those obtained in the control region. 

Interpolated threshold temperatures covering the surface of the extended region are 
then used in our algorithm for cloud detection in the B4 band. For a hxed instant, if the 
temperature obtained from the GOES image at 10.7 /im is colder than Tq then the presence 
of a cloud above the location is inferred. 


4. Threshold temperatures at the upper troposphere 

Clouds along the upper troposphere are detected following a similar schema, dehning 
hrst a threshold temperature for high altitudes where the weighting functions of B3 shows 
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the maximum value. We denote the altitude of the atmosphere measured from the ground 
with H, different to the altitude of a weather station h. We recall that the weighting 
functions depending on both, the wavelength and H, determine the layer from which the 
emitted radiation originates. For the B3 band, the median height values of the weighting 
functions are in the interval 8 — 10 km of altitude above the ground with a maximum at 9 
km (Cavazzani et al.2011). 

Threshold temperatures at upper troposphere may be obtained using different 
approaches. Hydayat et ah (2012) considered the reference atmospheric model of Kneizys 
et ah (1996) which provides a vertical tropospheric prohle of temperature for latitudes 
within ±15° and adopted a threshold equal to —30.15°C at the observatory location based 
on a analysis with data of sky visual inspection. Our approach is based on the use of 
three thresholds each one at heights of 8, 9 and 10 km obtained from radiosonde upper-air 
measurements. 


We make a strong emphasis on the reliability of the temperature at upper troposphere. 
Particularly, we are interested in hnding trustworthy threshold values between 8 and 10 
km. Data were obtained twice a day from the freely available database of the University 
of Wyoming which provides information of temperature, humidity, pressure and mixing 
ratio as a function of the altitude for two stations located inside the extended region, 80222 
SKBO (74°9’0”5W, 4°42’0”N, 2546m) in central Andes and 78988 TNCC (68°57’36”W, 
12°11’60”, 54m) in the Caribbean coast. Upper-air temperatures of both stations exhibit 
similar patterns as a function of the altitude. Around 60 radiosonde records were used per 
month adding data from both stations. Altitude resolution is variable with typical values 
ranging from 200 up to 600 m. 


^http://weather.uwyo.edu/upperair/sounding.html 
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We compute monthly mean temperatnres by considering measurements made twice 
daily (00 and 12 UTC) during the years 2008-12 at each H using linear interpolation along 
a uniform vertical grid of 21 points between 8 and 10 km. In order to follow the same 
approach than the one used for ground thresholds, we define the threshold temperature at 
the altitude H as Th =< Th > San where < Th > is the monthly mean value and an 
is the standard deviation. Monthly thresholds as a function of the altitude are displayed 
in Figure [6l Monthly standard deviations change with altitude and time showing yearly 
variations up to ±2°C in the worst of cases. In the tropics, the temperature anomalies at 
high altitudes is enhanced by the presence of moderate size clouds usually associated with 
shallow convection responsible for anomalies up to +2°C at altitudes from 5 up to 14km 
(Folkins 2008). This tropospheric warming has been largely observed and we suggest that 
is the main source of the sigma variations. 

The further assumption that these thresholds describe the upper troposphere of the 
whole extended region in many aspects is a strong hypothesis, however the area of study 
is affected by atmospheric phenomena including the westerly low level Choco jet coming 
from the Pacific coasts of Colombia (Poveda & Mesa 2000) and the Caribbean trade winds 
which leads to similar daily patterns of high altitude clouds. Middle and right panels in 
Figure [7| show a typical evolution during the night of 8th of April 2009. We see two images 
in B3 taken at the beginning and at the end of the night. High altitude clouds shown in 
green color in the left panel coming from the Pacihc ocean are blocked by the Andes system 
leaving clear zones in the Andes of Venezuela and over the Amazon zones. 

We note that the fact of extrapolating thresholds over all the extended region, the 
possible changes in the instrumentation and the differences between night and daytime 
observations may affect the observed temperatures and therefore introduce considerable 
noise in the computed thresholds. However, the monthly thresholds obtained with these 


15 


two stations are in agreement with the predictions of the tropical atmospheric model of 
Kneizys (1996) indicated with the solid line in FignreEl 


5. Methodology and night classification 

We implemented an algorithm for percentage of clear sky calcnlation dnring the night 
throngh a comparison pixel to pixel of GOES 12/13 interpolated temperatnres with monthly 
threshold valnes obtained from long-term records of weather stations and radiosonde 
databases. We constrncted fonr monthly threshold images, one for band B4 labeled as Tq 
and three for B3 labeled as Th each one at heights H = 8,9, and 10 km. Daily comparisons 
were done in both bands at 00:15:15, 06:45:13 and 09:45:14 UTC with similar approach to 
those described in Hidayat et al. (2012) and Erasmns & Sarazin (2002). First we nsed daily 
radiances in the B4 band to compnte brightness and actnal temperatnres. If the satellite 
derived temperatnre is colder than the monthly threshold at B4 the presence of middle high 
clonds above the pixel image at the instant i is inferred : 


if EB4,i > To —> dear 

otherwise —> covered (4) 

where Tq is a monthly constant depending only on the altitnde of the terrain. If a pixel 
is classified as covered in B4 it will be also covered in the B3 band. On the contrary, if 
the pixel is clear in B4 we still have to verify the presence of high altitnde cirrns clonds by 
doing a farther and similar analysis nsing the B3 band. The pixel is classihed as clear or 
covered in B3 at the altitnde H in the following way : 
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if < Th > —Tbs^i < Scth —y clear 

otherwise —> covered (5) 

where < > denotes the monthly mean temperature of the upper troposphere at 

height H and an the standard deviation at that height. We compute three values of clear 
sky fraction at the instant i and at each altitude H using inequalities (4) and (5) by 
counting the number of clears or covered pixels. The total clear sky fraction at altitude H 
is then obtained from the average of the three instants i along the night. One value of CSF 
is computed in this way at each altitude. 


5.1. Classification of the nights 

The nights in the extended region are classihed as spectroscopic or photometric at 
each altitude H based on the values of total clear sky fraction. We define the quality of the 
night, at the height H, as spectroscopic if the total clear sky fraction is larger than 33%. 
As we computed the total CSF with three values, spectroscopic nights are those with 1/3 
of the night free of clouds at least. Sites that fulhll this requirement are more suitable for 
spectroscopic observations than for photometric ones as long as spectra can be acquired 
quickly in comparison with photometeric data. We note that this dehnition is based only 
on the amount of clear sky and does not involve any sky quality indicator. Spectroscopic 
Night Simulations (hereafter, SNS) provide monthly averages of clear sky fraction at three 
different altitudes for each pixel of the image satisfying the condition that CSF is larger 
than 33%. 

Nights with photometric quality on the other hand, are those with total sky fraction 
larger than 67% corresponding to 2/3 of clear sky during the night at least. If this condition 
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is fulfilled in a particular site, we argue than at the altitude H will be possible to conduct 
photometric observations. Photometric Night Simulations (hereafter, PNS) provide monthly 
maps of clear sky fraction at three different altitudes for each pixel of the image with the 
condition that CSF is larger than 67%. 


6. Validation of the algorithm 

Based on total CSF values obtained with SNS and PNS outputs, monthly averages at 
altitude H were computed. Cloud cover along the upper-troposphere is sensitive to the 
altitude, being larger at low altitudes mostly at H = 8 km, indicating that high cirrus 
does not affect considerably the net cloud fraction. We used the monthly CSF averages at 
H = 8, 9 and 10 km to compute a weigthed mean for the three altitudes at each instant 
during the night. Weights were computed through a comparison of the amount of clear 
sky fraction for photometric nights coming from our model with those obtained from daily 
number of observing hours reported by log-books of CID or Site 2 in Venezuela. Site 2 is 
placed at 3600 meters height and enclosed by higher mountains as Pico Bolivar (4978m) 
at less than 40 km away. As the resolution is 28 km in longitude and 18 km in latitude, 
the threshold temperature at Site 2 is highly affected by the altitude of the closest higher 
mountains giving colder thresholds than expected. Using PNS instead of SNS makes more 
precise validation because PNS output is closer to the condition of total clear night. 

Figure E] corresponds to the clear hours above the reference Site 2 or CID as reported 
in the log-book. Despite the bias present in the log-books and the missing data for the 
years 2009-10, due mainly to bad weather, we did not hnd considerable differences between 
observed and expected monthly clear sky fraction values when different weights were used. 
Therefore, for the purposes of this study, we consider sufficient to use equal weights for the 
three altitudes. 
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The validation of our code is limited due to factors such as the lack of completeness 
of daily log-books, tecnhical issues, high humidity, strong winds and uncertainties of our 
model. We used all the information available in the log books as the general description 
of the night done by the astronomer to compute the amount of daily clear sky fraction for 
each month of the years 2008-12 at Site 2. Comparison with our monthly results using PNS 
is displayed in Figure O The solid line represents no differences between predicted and 
measured monthly clear sky fraction values. There is a large spread during the considered 
years, especially for clear sky fraction values between 20-60%, but in general terms they are 
in reasonable agreement. 


7. Monthly clear sky fraction for 2008-12 

Monthly percentages of clear sky for the years 2008-12 over the extended region were 
obtained using the method described above. The main common feature is that monthly 
averages display maxima at the same locations along the hve years. Maximum monthly 
percentages of clear sky are observed during two periods in the years 2008-12, one from 
December to February and another from June to July according with the bimodal annual 
cycle of precipitation (Poveda et ah 2011). The year 2012 stands out as the best year 
displaying the highest percentages during both periods. 

We take the year 2012 as representative of the monthly distributions of clear sky. As 
illustrated in Figure fTOl our monthly averages of CSF for SNS nights exhibit differences 
along the year with the highest CSF values always located over a wide region in the central 
Venezuelan Andes with a maximum CSF>80% at Site 1 during both periods. For Site 2 or 
CID corresponding to Observatorio Nacional de Llano del Hato we get a CSF mean value 
of 58.5%, lower in comparison with Site 1 during both periods. 
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We have found two regions in Colombia with similar night qualities to those observed 
at Site 2. The northern part of the Sierra Nevada de Santamarta close to Nabusimake 
(NAB) and the central Serrania del Perijd toward the east of Codazzi (COD). Percentages 
at these two regions are in the range of 80-90% during the first period, 50-70% during the 
second and below 30% from March to May and between August and October. Sites over 
this region contain the best places in Colombia found in this study. 

Somewhat low quality places are found over the north-eastern Andes of Colombia, over 
mountain systems along BUC, MAC and PIS with clear sky fraction between 30-50% during 
both periods of 2012. Southern Andes regions in the Nudo de los Pastes, especially those 
in northern Ecuador display CSF values below 30% from November to May and ~70% 
from June to October. Very low percentages of clear sky fraction below 20% characterize 
central Andes of Colombia over the regions comprising Popaydn (POP) and BOG with the 
exception of the region centered at the high site Nevada del Huila (75°55’12”W, 2°51’36”N) 
at 3700 meters with a mean of 35% between December-February and June-July and 33% 
for the other months. 

Regarding monthly averages of CSF for PNS nights, the same regions were found but 
with low percentages in comparison as shown in Figure [TTJ The absence of spectroscopic 
nights in central Colombia from May to October conhrms that November-February and 
June-July are marked by the largest CSF percentages while March-May and August- 
November are characterized by lower CSF values. We also observe such behaviour for other 
years. The existence of this bimodal behaviour with maxima during the dry season and 
minima during the rainy season is atributed, as we metioned already, to the double passage 
of the intertropical convergence zone for the extended region. Therefore we conducted yearly 
SNS and PNS simulations, bringing together the months December, January, February, 
June, July and August into one group that we call dry season and March, April, May, 
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September, October and November into another gronp denoted as rainy season. 


SNS percentages for both seasons show maxima at the same locations than those 
observed for the year 2012 with seasonal differences of ~20% with the exception of low 
latitnde locations close to the eqnator where slightly smaller differences are observed. 
Percentages for PNS simnlations resnlt 10% lower in comparison to SNS nights bnt with 
identical CSF distribntion over the extended region. 


Annnal clear sky fraction averages obtained from onr SNS simnlations dnring rainy 
seasons for the years 2009 and 2012 are larger than for the other years. It seems to correlate 
with the changes in the Oceanic Nino Index (ONI) reported by the climate prediction center 
NOAAo. Changes from positive to negative anomalies in April 2010 associated to a La Nina 
event increase the nnmber and intensity of tropical storms in central Colombian Andes thns 
decreasing the clear sky percentages. The Nudo de los Pastos and the Sierra Nevada de 
Santamarta are not considerably affected by this El Nino Southern Oscillation (ENSO). 


8. Selected sites from five year CSF averages 

For selection pnrposes we give preference to the simnlations for spectroscopic nights 
dnring the dry seasons. We have selected twelve sites based on Eve year averages dnring 
the dry seasons that can be snitable for astronomical observations by selecting locations 
with CSF>30%. Factors snch as accesibility, secnrity, service, local hnmidity, atomospheric 
seeing were not considered and we need additional stndies to conhrm these locations as 
potential for astronomical sites. These twelve selected sites are distribnted mostly along 
the following fonr zones: the central Andes of Veneznela comprising the states of Merida 
and Trujillo, the north east of the Sierra Nevada de Santamarta aronnd sites NAB and 


“^http://www. cpc.noaa.gov/prodncts/analysis_monitoring/ensostnff/ensoyears.shtml 
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COD, the north eastern Andes of Colombia in the proximity of El Espino (ESP) and a wide 
region on the southern Andes in Nudo de los Pastos containing two locations in Ecuador. 
We show in Figure [T^ the contours of clear sky fraction for these selected zones during the 
dry season with a spatial pixel resolution of 28 x 18 km. 

Even though the selection criteria involves only spectroscopic nights during the dry 
seasons, we give a total number of clear nights per year at selected sites by adding the CSF 
hve year averages for each season and then multiplying by a factor of 1.827. Tables [2] and 
|3] contain properties of the selected sites indicating the identifier, the nearest populated 
location, geographic coordinates of the site, altitude, yearly hve year averages for the dry 
and rainy seasons, hve year average of clear sky fraction, standard deviation and the total 
number of clear nights per year. 


8.1. Central Venezuelan Andes 

The best places are situated over an extense area of the central Andes of Venezuela 
comprising latitudes from 8°48’N to 9°11’N and longitudes between 70°17’W and 70°42’W. 
The region contains a number of potential sites at high altitudes, in particular Site 1 
(70°28’48”W, 9°5’60”N) where we obtain the maximum CSF in the extended region. 

Site 1 has an altitude of 3480 meters and it is located ~35 km north-east from Timotes 
(TIM) very close to the protected area Parque Nacional General Cruz Carrillo. Five year 
average of clear sky fraction for spectroscopic nights during the dry seasons for this site 
is (67.0±10.7)% corresponding to 122.3 clear nights per year. The percentage during the 
rainy seasons is (53.5±12.2)% leading to 97.7 clear nights per year. Subsequently we report 
a total number of (220±30) spectroscopic nights per year at Site 1. 


Regarding our reference Site 2 or CID, hve year average of spectroscopic nights during 


22 


dry seasons is (46.9±8.1)% while for rainy seasons it is (25.5±10.8)% giving a mean of 72.4% 
corresponding to (132±25) total clear nights per year. Five year average for photometric 
nights at this site gives (28.9±3.4)% for dry seasons and (14.1±8.6)% for rainy seasons 
giving a total of 43% which is in agreement with the log-books of the observatory as shown 
in Fignre [9] even thongh seasonal differences were not considered for the validation of the 
algorithm as discnssed in Section 6. 


8.2. Sierra Nevada de Santamarta 

This is a large system of monntains in the north of Colombia isolated from the Andes 
and located jnst 52 km from the Caribbean coast. It rises np to 5700 meters above sea 
level over a protected area of 17000 km^ containing the national park Sierra Nevada de 
Santamarta. Five year average of CSF is (42±10)% over a broad region of 40 km to the 
north of NAB, specifically at Site 3 (73°31’12”W, 10°55’12”N) placed at an altitnde of 3600 
meters with 76.8 clear nights per year dnring the dry season and 49.3 nights per year dnring 
the rainy season. The total spectroscopic clear nights at Site 3 is (126±34) in contrast with 
(72±16) clear nights per year with photometric qnality. 

The steep gradient in elevation toward the sonth-east direction of NAB make the 
percentage of CSF lower along the plain regions on the valley of river Cesar which separates 
the Sierra Nevada de Santamarta from the Serrama del Perijd, an extension of the 
eastern Andean branch on the bonndary between Colombia and Veneznela. Slightly lower 
percentages than those observed at Site 3 are obtained at Site 4 less than 10 km from the 
location COD. The site has (111±27) clear nights per year with spectroscopic qnality and 
(66±18) photometric nights. 
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8.3. Eastern Andes of Colombia 

This broad territory hosts a high number of Paramos, dry forests, valleys and high 
snow peaks most ot them protected areas such as national parks El Cocuy (COC) and 
Pisba (PIS). The selected locations are within the wide region between the Sierra Nevada 
del Cocuy and the north part of the canyon of the river Chicamocha. The highlands of the 
sierra reach an elevation of more than 5000 meters including several snow peaks. Site 6 
placed ~50 km west from these high snow peaks, is the center of a region with five year 
average for the dry seasons of (35±4)%. Clear sky percentage toward the north west along 
the canyon reaches the maximum value of (37.4±5.7)% at Site 5 located in Cuchilla las 
Torrecitas, next to Cordillera de Hoya Crande in Norte de Santander at an altitude of 4060 
meters. For rainy seasons we obtain (21.3±6.7)% at this site giving a total of (107±16) 
clear nights per year. 

Site 7 located 20 km north-west of PIS with an altitude of 3600 meters belongs to 
the protected area Parque Nacional Paramo de Pisba. Five year average is (34.8±2.9)% 
whereas for rainy seasons Site 7 has 25% giving a total of (99±16) clear nights per year. 
Another selected site in this region is Site 8, the closest selected location to the snow peaks, 
with (90±12) clear nights per year using only four years since the 2012 CSF average is 
below our detection threshold. Percentages of clear sky fraction for photometric nights for 
Sites 7 and 8 are indicated in Table [3l 

Lower quality sites with CSF<15% are clearly observed over the central branch of the 
Andes comprising locations Murillo (MUR) and Medellin (MED). Brieva (1985), using 
cloud cover data from visual inspection of the sky reports between 130 and 250 useful nights 
per year at site Otun (4°46’33”N, 75°24’53”W), located 25 km in the south west direction 
from MUR. We obtain (89±11) clear nights per year at site MUR (4°59T9”,75°10T2”). 
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8.4. Nudo de los Pastas 

This region consists of tropical hnmid forests where the Andes splits into the western 
and central branches. Sites 9 and 11 placed in north of Ecuador stand out as the best 
places in the region with (105± 36) and (101± 42) clear nights per year, respectively. With 
(24.3±3.9)% below our selection limit we include the Site 12 Laguna de la Cocha (LCO) 
located at the Colombian side of the Nudo de los Pastos. 

Minimum seasonal differences are found at the sites in this region. Comparatively with 
the other sites, the percentage of clear sky during dry and rainy seasons is very similar 
in sites within the Nudo de los Pastos confirming that the sky cover does not change 
signihcantly from one season to another at locations near the equator. We note from Tables 
[2] and [3] that the minimum differences between CSF values for dry and rainy seasons for all 
twelve sites are observed especially for the sites 9, 11 and 12 with differences of 4, 2 and 1%. 

Five year average for the selected sites exhibit an extended minima in 2010-11 during 
rainy seasons as shown in Figure [131 We found a possible correlation to La Nina episodes 
since positive values for the Southern Oscillation Index (SOI) are observed for these two 
years. Interior plots in Figure IT^ show SOI values obtained from NOAA website for 2008-12. 
The symbols for the months go along with those for the sites. Positive values mostly for 
2008, 2010 and 2011 are associated to intensification of storms over the Andes and then to 
low CSF values as observed. 


9. Conclusions 

In this work we describe a method to compute the clear sky fraction from GOES data 
based on comparisons with long-term records of temperature from weather stations and 
radiosondes, over a wide region of northern South-America comprising sites in Colombia, 


Venezuela and Ecuador. 


We conducted a temporal analysis based on three daily values along five years doing 
a careful comparison of GOES12/13 temperatures in bands B3 and B4 with thresholds at 
different altitudes. This approach enabled us to carefully determine the percentage of clear 
sky simultaneosly over the whole region. Similar perentages were obtained for all locations 
along the five years confirming the reliability of our methodology. 

Threshold monthly temperatures as a function of the altitude of the terrain were 
obtained for a small region of control through Kriging interpolation of long-term records 
of temperature at weather stations. Monthly linear correlations between thresholds and 
altitude were found and extrapolated to an extended region including Observatorio Nacional 
de Llano del Hato. Threshold monthly temperatures at 8, 9 and 10 km along the upper 
atmosphere were also computed based on daily radiosonde data of two stations located at 
sea level and at the central Andes. 

The method differs from other cloud fraction algorithms that uses long-term records of 
temperature to obtain the CSF. Our approach gives more precise CSF values as long as 
we have used reliable reference temperatures at the surface and at the high troposphere. 
Furthermore, our model can be easily extended to other tropical regions of astronomical 
interest. 

Nights are classified as spectroscopic or photometric based on the three CSF values 
during the night. The total number of spectroscopic and photometric nights per year at a 
fixed site are computed separately from five year averages for dry and rainy seasons. 

The best location found in this study is Site 1 (70°28’48”W, 9°5’60”N) in the central 
Venezuelan Andes with (220±42) spectroscopic clear nights per year and (157±24) 
photometric clear nights per year. Site 1 is placed at an altitude of 3480 meters and 35 km 
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north-east from Timotes within a protected area of Venezuela. Two lower quality locations 
were found in Sierra Nevada de Santamarta, Site 3 (73°31T2”W, 10°55T2”N) and Site 4 
(73°2’24”W, 10°8’24”N) with (126±34) and (111±27) spectroscopic clear nights per year, 
respectively. Site 3 stands out as the best site identihed in Colombia. 

The eastern Andes of Colombia host two regions with more than 90 spectroscopic 
clear nights per year. The region centered at Site 5 (72°43T2”W, 7°1T2”N) in Norte de 
Santander with (107±23) clear nights per year and the wide zone centered at El Espino 
comprising Macaravita and El Cocuy in Boyacd with (94±13) clear nights per year. The 
nearest locations to the equator in the Nudo de los Pastos region are Site 9 ( 77°50’24”W, 
0°46’48”N) and Site 11 (77°50’24”W, 0°15’36”N) with (105±36) and (101±42) clear nights 
per year, respectively, presenting minimum seasonal differences of number of clear nights as 
expected. Finally, CSF values below 25% observed during the rainy seasons of the years 
2010 and 2011 seems to be correlated to positive values of the Southern Oscillation Index. 

The selected sites obtained in this study may be considered as a hrst-step to decide 
the location of an astronomical observatory. However, site-survey is absolutely required 
before start any astronomical facility. Results of an observational campaign with the aim 
to characterize the sites identihed in this study will be reported in a forthcoming paper. 

We would like to further point out that the value of our approach goes beyond the 
possibility of being useful for identihcation of potential sites in the region. Our approach 
allows us to quantify the contribution of the different layers of the atmosphere to the CSF 
above the site, in particular, the high troposphere and the surface. Although we have had 
to make a number of simplifying assumptions, the results are self-consistently and giv 
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Fig. 2.— Historical monthly thresholds of temperature at the locations of the weather 
stations with long-term records as a function of the altitude. Solid line represents the trend 
of the minima temperature obtained from short-term records: Tmin,tr = —5.5°C/kmxh -|- 
23.42°C where h is the altitude of the station in km. 



Semivariogram 
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Di stance 


Fig. 3.— Semivariograms computed along south-north (hlled circles) and west-east directions 
(open circles). Both semivariograms are well described by a spherical model (dashed line) 
with a nugget of 0.048, sill of 1.01 and range of 120. 
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Fig. 4.— Monthly threshold temperatures for January in the control region. Panel a) 
Thresholds obtained through a Delanuay triangulation of the IDEAM long-term records 
(open squares), b) thresholds using Kriging interpolation of 400 altitudes from IDEAM all¬ 
term records, c) sigma errors of the Kriging interpolation and d) cross validation of Kriging 
altitudes with IDEAM values. A perfect interpolation is indicated with a solid line. 
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Fig. 5.— Monthly threshold temperatures for cloud detection in the B4 band over the ex¬ 
tended region for April and December which are representative of the rainy and dry seasons, 
respectively. The highest mountains in branches of the Andes shown in blue present thresh¬ 
olds below 5°C. Caribbean and Pacihc coasts, the valleys of the Magdalena river and the 
flat south-eastern region are marked by thresholds between (20-25)°C. 
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Month 

A (°C/km) 

B (°C) 

JAN 

-5.98±0.12 

24.48±0.27 

FEB 

-6.01±0.13 

24.42±0.29 

MAR 

-5.96±0.12 

24.71T0.27 

APR 

-5.83±0.11 

24.85±0.25 

MAY 

-5.71±0.11 

24.60T0.26 

JUN 

-5.69±0.13 

24.34±0.28 

JUL 

-5.93±0.15 

24.47±0.33 

AUG 

-5.84±0.14 

24.41±0.30 

SEP 

-5.79±0.12 

24.35±0.26 

OCT 

-5.72±0.11 

24.35±0.24 

NOV 

-5.66±0.11 

24.28±0.26 

DEC 

-5.79±0.13 

24.22±0.28 


Table 1: Historical monthly lapse rates (H) and intercepts {B) obtained from linear regres¬ 
sion of threshold temperatures (from Kriging) with the altitude of the weather station (i.e. 
To(h) = A X h + B). Variations of 0.1°C/km in the lapse rate are significant for threshold 
computations whereas percentage errors in B are of the order of 1% and are considered 
insignihcants. 
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Fig. 6.— Monthly threshold temperatures for cloud detection in the B3 band over the 
extended region for March and April which are representative of the rainy season and for 
and December and January representative of the dry season. Yearly standard deviations in 
the upper atmosphere obtained from radiosonde data are also displayed in the interior plot. 
Temperature from the tropical atmospheric model of Kneizys et ah (1996) is indicated with 
a solid line. 
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Fig. 7.— Temperature at upper atmosphere using B3 radiances. High altitude clouds in 
blue move from left to right in timescales of days before being blocked by the Andes system, 
leaving clear regions mostly in the western Venezuela and in the eastern Andes of Colombia. 
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Fig. 8.— Number of daily observing hours obtained from the log-books of Observatorio 
Nacional de Llano del Hato (CID or Site 2) for the years 2008-12. 
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Fig. 9.— Absolute difference between clear sky fraction from our model and those obtained 
from long-term records of the number of observing hours for 2008-12 at CID or Site 2. A 
linear regression leads to a slope of (0.76±0.13) corresponding to a percentage difference of 
24% in comparison with the perfect agreement between observed and expected CSF values. 
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Fig. 10.— Monthly clear sky fraction percentage for spectroscopic nights (SNS) in the 


extended region for the year 2012 for April and December which are representative of the 
rainy and dry seasons, respectively. 
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Fig. 11.— Monthly clear sky fraction percentage for photometric nights (PNS) in the ex¬ 


tended region for the year 2012 for April and December which are representative of the rainy 
and dry seasons, respectively. 
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Fig. 12.— Five year average of clear sky fraction for spectroscopic nights during the dry 
seasons in four regions containing the twelve selected sites in this study. 
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Fig. 13.— Yearly clear sky fraction for spectroscopic nights for twelve selected sites for 
dry (left) and rainy (right) seasons of the years 2008-12. Dry Season comprises the months 
of December, January, February, June, July and August and Rainy Season the months of 
March, April, May, September, October and November. Interior plot corresponds to the 
Southern Oscillation Index (SOI) in the Nino 3.4 region (-5S,5N;-120-170W) for the same 
years. Symbols of the monthly anomalies go along with the numbered sites. 







Table 2. Percentages of clear sky fraction for spectroscopic nights (PNS) for selected sites. Dry Season comprises the 
months of December, January, February, June, July and August and Rainy Season the months of March, April, May, 

September, October and November. 


Site 

ID 

Nearest 

Location 

Lon. 

(W) 

Lat. 

(N) 

Alt. 

(m) 

2008 

Dry Rainy 

2009 

Dry Rainy 

2010 

Dry Rainy 

2011 

Dry Rainy 

2012 

Dry Rainy 

Mean CSF(%) 

Dry Rainy 

No.nights 

per year 

1 

TIM 

Timotes 

70°28’48” 

9°5’60” 

3480 

65.4 

52.1 

68.8 

56.8 

55.5 

42.7 

61.1 

43.4 

84.0 

72.5 

67.0±10.7 

53.5±12.2 

220±42 

2 

CID 

Cida 

70°52T2” 

8047-24” 

3600 

51.4 

27.0 

45.0 

26.6 

38.7 

20.2 

40.8 

12.3 

58.5 

41.6 

46.9±8.1 

25.5±10.8 

132±34 

3 

NAB 

Nabusimake 

73°31T2” 

10°55T2” 

3300 

43.4 

29.8 

40.4 

36.2 

33.0 

17.8 

35.3 

17.9 

58.1 

33.3 

42.0±9.9 

27.0±8.7 

126±34 

4 

COD 

Codazzi 

73°2’24” 

10°8’24” 

2620 

42.6 

26.7 

38.0 

31.8 

33.5 

17.0 

28.6 

13.1 

46.7 

26.2 

37.9±7.2 

23.0±7.7 

111±27 

5 

CHI 

Chitaga 

72°43T2” 

701-12- 

4060 

39.5 

20.8 

36.6 

25.7 

31.7 

16.0 

33.1 

14.0 

45.9 

30.2 

37.4±5.7 

21.3±6.J 

107±23 

6 

ESP 

El Espino 

72°24’0” 

6°30’0” 

3620 

33.6 

19.9 

37.5 

20.3 

29.5 

14.6 

34.4 

14.0 

39.9 

13.5 

35.0±4.0 

16.5±3^ 

94±13 

7 

PIS 

P. Pisba 

72°33’36” 

5°58’48” 

3600 

37.4 

19.6 

35.4 

25.4 

32.0 

15.5 

31.7 

12.2 

37.7 

25.5 

34.8±2.9 

19.6±5.^ 

99±16 

8 

COC 

El Cocuy 

72°4’48” 

6°30’0” 

3600 

31.8 

17.2 

35.2 

21.7 

29.8 

12.9 

35.5 

14.0 



33.1±2.8 

16.5±3.9 

90±12 

9 

TUF 

P. Tuffino 

77°50’24” 

0°46’48” 

3400 

26.6 

12.2 

31.8 

39.5 

27.9 

15.7 

28.0 

22.0 

39.6 

44.5 

30.8±5.3 

26.8±14.4 

105±36 

10 

HUI 

N. Huila 

75°55T2” 

2°51’36” 

3700 

27.6 

12.4 

28.2 

24.8 

34.0 

10.5 

24.9 

17.8 

35.4 

33.2 

30.0±4.5 

19.7±9.4 

91±25 

11 

IBA 

Ibarra 

77°50’24” 

0°15’36” 

3660 

20.3 

10.8 

33.4 

42.5 

22.1 

14.5 

26.9 

23.8 

41.0 

40.7 

28.7±8.5 

26.5±14.6 

101±42 

12 

LCO 

L. Cocha 

77°2’24” 

102-24” 

3400 

18.3 

9.3 

23.3 

32.01 

28.2 

13.1 

24.5 

17.3 

27.0 

42.1 

24.3±3.9 

22.8±13.8 

86±32 






Table 3. Percentages of clear sky fraction for photometric nights for selected sites. Dry Season comprises the months 
of December, January, February, June, July and August and Rainy Season the months of March, April, May, 

September, October and November. 


Site 

ID 

Nearest 

Location 

Lon. 

(W) 

Lat. 

(N) 

Alt. 

(m) 

2008 

Dry Rainy 

2009 

Dry Rainy 

2010 

Dry Rainy 

2011 

Dry Rainy 

2012 

Dry Rainy 

Mean CSF(%) 

Dry Rainy 

No.nights 

per year 

1 

TIM 

Timotes 

70°28’48” 

9°5’60” 

3480 

54.3 

36.7 

55.5 

35.2 

44.2 

29.1 

44.8 

28.5 

57.9 

44.8 

51.3±6.4 

34.9±6.6 

157±24 

2 

CID 

Cida 

70°52T2” 

8°47’24” 

3600 

29.6 

14.1 

29.3 

13.1 

25.8 

13.1 

25.7 

3.1 

34.0 

27.3 

28.9±3.4 

14.1±8.6 

79±22 

3 

NAB 

Nabusimake 

73°31T2” 

10°55T2” 

3300 

31.8 

17.9 

24.2 

16.9 

24.3 

8.7 

24.2 

7.1 

30.5 

11.4 

27.0±3.8 

12.4±4.8 

72±16 

4 

COD 

Codazzi 

73°2’24” 

10°8’24” 

2620 

29.5 

15.8 

26.7 

13.8 

22.3 

9.3 

16.5 

4.9 

28.7 

12.3 

24.7±5.4 

11.2±4.3 

66±18 

5 

CHI 

Chitaga 

72°43T2” 

7°1T2” 

4060 

27.6 

10.8 

23.2 

11.4 

20.7 

8.2 

20.7 

6.7 

29.3 

10.9 

24.3±4 

9.6±2.0 

62±11 

6 

ESP 

El Espino 

72°24’0” 

6°30’0” 

3620 

22.0 

7.6 

23.4 

8.1 

23.2 

8.8 

24.0 

8.1 

26.2 

4.8 

23.8±1.5 

7.5±1.6>F 

57±6 

7 

PIS 

P. Pisba 

72°33’36” 

5°58’48” 

3600 

25.4 

10.3 

23.2 

12.2 

23.2 

6.3 

21.3 

5.6 

26.8 

6.4 

24.0±2.1 

8.2±2.9 

59±9 

8 

COC 

El Cocuy 

72°4’48” 

6°30’0” 

3600 

21.6 

6.4 

23.1 

9.9 

21.7 

7.2 

22.1 

7.3 



22.1±0.7 

7.7±1.5 .1 

54±4 

9 

TUF 

P. Tuffino 

77°50’24” 

0°46’48” 

3400 

6.8 

6.5 

18.6 

17.5 

15.6 

8.5 

15.2 

9.6 

25.1 

28.6 

16.3±6.6 

14.1±9.1 

56±29 

10 

HUI 

N. Huila 

75°55T2” 

2°51’36” 

3700 

12.3 

5.1 

10.7 

10.6 

19.9 

6.0 

13.3 

8.7 

13.5 

8.8 

13.9±3.5 

7.8±2.2 

40±11 

11 

IBA 

Ibarra 

77°50’24” 

0°15’36” 

3660 

7.9 

6.0 

21.1 

17.8 

10.3 

8.0 

15.5 

14.6 

22.7 

29.0 

15.5±6.5 

15.1±9.1 

56±29 

12 

LCO 

L. Cocha 

77°2’24” 

1°2’24” 

3400 

5.3 

2.7 

11.0 

12.1 

14.3 

4.0 

11.3 

7.3 

13.3 

22.8 

11.0±3.5 

9.8±8.1 

38±21 






